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Abstract. We investigate the problem of reheating and its effects on the evolution of 
the curvature perturbations in tachyonic inflationary models. We derive the equations 
governing the evolution of the scalar perturbations for a system consisting of a tachyon 
and a perfect fluid. Assuming the perfect fluid to be radiation, we solve the coupled 
equations for the system numerically and study the evolution of the perturbations 
from the sub-Hubble to the super-Hubble scales. In particular, we analyze the effects 
of the transition from tachyon driven inflation to the radiation dominated epoch on 
the evolution of the large scale curvature and non-adiabatic pressure perturbations. 
We consider two different potentials to describe the tachyon and study the effects of 
three possible types of decay of the tachyon into radiation. We plot the spectrum 
of curvature perturbations at the end of inflation as well as at the early stages of 
the radiation dominated epoch. We find that reheating modifies the amplitude of 
the curvature perturbations only when the decay rate of the tachyon into radiation 
is explicitly time-dependent. We illustrate that it is the growth of the relative non- 
adiabatic pressure perturbation between the tachyon and radiation that is responsible 
for the change in the amplitude of the curvature perturbations during the transition. 
Our results corroborate similar conclusions that have been arrived at earlier based 
on the study of the evolution of the perturbations in the super-Hubble limit. In an 
appendix, we also discuss the corresponding results for the popular chaotic inflation 
model in the case of the canonical scalar field. We briefly comment on the possible 
implications of our results. 
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1. Effects of inflationary and post-inflationary dynamics on the evolution of 
the super-Hubble perturbations 

In the standard inflationary scenario, the scales of cosmological interest exit the 
Hubble radius within the first few (about 8-10) e-folds of inflation and, hence, are 
outside the Hubble scale during the later epochs. When comparing the predictions 
of the inflationary models with the Cosmic Microwave Background (CMB) and the 
large scale structure data, it is often assumed that the amplitude of the curvature 
perturbations remain constant at super-Hubble scales. In such a situation, the observed 
CMB anisotropies directly determine the amplitude and the shape of the perturbation 
spectrum imprinted on the modes when they left the Hubble radius during inflation. 
Typically, the amplitude of the perturbations constrain the parameters that describe 
the inflaton potential, while the shape of the spectrum limits its form (see any of the 
standard texts [1] or one of the following reviews [2]). 

Provided inflation is of the slow roll type for all of the required number of e-folds, 
it is indeed true that the amplitude of the curvature perturbations freeze at their value 
at Hubble exit. However, if there is a period of deviation from slow roll inflation, then 
the asymptotic (i.e. the extreme super-Hubble) amplitude of the modes that leave the 
Hubble radius just before the deviation are enhanced when compared to their value at 
Hubble exit [3]. While modes that leave well before the deviation remain unaffected, it 
is found that there exists an intermediate range of modes whose amplitudes are actually 
suppressed at super-Hubble scales [4, 5]. Depending on the form of the departure from 
slow roll, these effects lead to certain features in the scalar perturbation spectrum. 
If such deviations occur either during or soon after the cosmological scales leave the 
Hubble radius, then the CMB observations constrain the resulting features rather well 
(for an inexhaustive list, see Refs. [6]). But, at smaller scales, only theoretical tools are 
currently available to restrict the form of the primordial spectrum. These constraints 
are essentially based on the number density of primordial black holes that are formed 
towards the end of inflation (for recent discussions in this context, see, for example, 
Refs. [7]) 

Since the cosmological scales are well outside the Hubble radius by the early stages 
of inflation, clearly, the shape of the perturbation spectrum on such large scales is indeed 
unlikely to be affected by subsequent dynamics. But, over the past decade, it has been 
recognized that post-inflationary dynamics can alter the amplitude of the curvature 
perturbations at super-Hubble scales, in particular, when more than one component 
of matter is present. Preheating, the curvaton scenario and the modulated reheating 
mechanism are popular examples that illustrate the interesting possibilities of post- 
inflationary dynamics. Preheating — a mechanism that transfers energy from the inflaton 
to radiation through an explosive production of quanta corresponding to an intermediate 
scalar field — is known to even lead to an exponential growth in the amplitude of the 
super-Hubble perturbations (for the earlier discussions, see Refs. [8] and, for more recent 
efforts, see, for instance, Refs. [9]). In the curvaton scenario, while the inflationary epoch 
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still remains the source of the perturbations, these perturbations are amplified after 
inflation due to the presence of entropic perturbations (see, for example, Refs. [10]). The 
modulated reheating scenario is an extreme case wherein inflation is essentially required 
only to resolve the horizon problem, whereas the perturbations are generated due to 
an inhomogeneous decay rate when the energy is being transferred from the inflaton 
to radiation through other fields [11, 12, 13]. These different alternatives indicate that 
the post-inflationary evolution of the large scale curvature perturbations can be highly 
model dependent and, therefore, requires a careful and systematic study. Evidently, in 
these scenarios, the effects of post-inflationary dynamics have to be taken into account 
when constraining the inflationary models using the CMB observations (in this context, 
see, for instance, Ref. [14]). 

With these motivations in mind, in this paper, we investigate the problem of the 
more conventional reheating scenario [15, 16, 17], and its effects on the evolution of the 
curvature perturbations in tachyonic inflationary models (for the original discussions 
on the tachyon, see Refs. [18, 19]; for efforts on treating the tachyon as an inflaton, 
see, for instance, Refs. [20, 21, 22]; for discussions on reheating in such inflationary 
models, see Refs. [23, 24]). We shall consider two types of potentials to describe the 
tachyon and construct scenarios of transition from inflation to radiation domination for 
the following three possible cases of the decay rate T of the tachyon into radiation: (i) a 
constant, (ii) dependent on the tachyon, and (iii) explicitly dependent on time. We 
solve the coupled equations for the system numerically and study the evolution of the 
perturbations from the sub-Hubble to the super-Hubble scales. Importantly, we shall 
consider the effects of the transition from inflation to the radiation dominated epoch 
on the evolution of the large scale (i.e. those that correspond to cosmological scales 
today) curvature and non-adiabatic pressure (i.e. the intrinsic entropy as well as the 
relative entropy or the isocurvature) perturbations. We shall evaluate the spectrum 
of curvature perturbations at the end of inflation as well as at the early stages of the 
radiation dominated epoch. As we shall illustrate, reheating affects the amplitude of 
the curvature perturbations only when the decay rate of the tachyon into radiation is 
explicitly time dependent. We shall also show that it is the growth of the isocurvature 
perturbation (i.e. the relative non-adiabatic pressure perturbation between the tachyon 
and radiation) that induces the change in the amplitude of the curvature perturbations 
during the transition. Our results confirm similar conclusions that have been arrived at 
earlier in the literature based on the study of the evolution of the perturbations in the 
super-Hubble limit [11]. We shall also discuss similar effects in the case of the canonical 
scalar field. 

A few clarifying remarks are in order at this stage of the discussion. We should stress 
that the set up we are considering is the standard cold inflationary scenario, followed 
by an epoch of reheating achieved by the standard method of introducing a coarse- 
grained decay rate in the equation of motion describing the inflaton [15, 16]. Recently, 
an analysis somewhat similar to what we shall consider here has been studied in the 
context of warm inflationary scenarios involving tachyonic fields [25]. Though there can 
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exist similarities in the form of the equations in the cold and the warm inflationary 
scenarios, the physics in these two scenarios are rather different (for a recent discussion, 
see, for example, Ref. [26]). Moreover, to analyze the effects of reheating on the large 
scale perturbations, one could have possibly studied the evolution of the perturbations 
in the super-Hubble limit, say, the first order (in time) differential equations usually 
considered in the literature (cf. Refs. [11, 12, 27]). However, some concerns have been 
raised that the coupling between the inflaton and radiation may affect the amplitude 
as well as the scalar spectral index in certain situations [16, 17]. Also, the modified 
background dynamics preceding reheating can effect the extent of small scale primordial 
black holes that are formed towards the end of inflation. If we are to address such issues, 
it requires that we study the evolution of the perturbations from the sub-Hubble to the 
super-Hubble scales. We shall comment further on some of these points in the concluding 
section. 

The remainder of this paper is organized as follows. In the following section, 
we shall summarize the essential background equations for the system consisting of 
the tachyon and a perfect fluid and set up scenarios of transition from inflation to 
radiation domination. In section 3, we shall obtain the equations describing the scalar 
perturbations for the system of the tachyon and a perfect fluid. In section 4, we shall 
evolve the coupled system of equations describing the perturbations and study the effects 
of the transition from inflation to the radiation dominated epoch on the super-Hubble 
curvature perturbations. We shall plot the spectrum of curvature perturbations at the 
end of inflation as well as at the early stages of the radiation dominated epoch. We 
shall also explicitly illustrate that it is the growth of the isocurvature perturbations 
that acts as the source of change in the amplitude of the curvature perturbations during 
the transition. Finally, in section 5, we conclude with a summary and a discussion on 
the results we have obtained. In the appendix, we shall discuss the corresponding results 
for the chaotic inflation model in the case of the canonical scalar field. 

The notations and conventions we shall adopt are as follows. We shall set h — c — 1, 
but shall display G explicitly, and define the Planck mass to be M P1 = (8ttG)" 1/z . We 
shall work with the metric signature of (+,—,—,—). We shall express the various 
quantities in terms of the cosmic time t, and we shall denote differentiation with respect 
to t by an overdot. While we shall use the subscript i to denote the spatial components 
of, say, the momentum flux associated with the matter fields, the subscripts a and (3 
shall refer to the two components of matter. 

2. The transition from tachyon driven inflation to the radiation 
dominated epoch 

In this section, we shall discuss the background equations describing the evolution of 
the tachyon that is interacting with a perfect fluid. Assuming the perfect fluid to 
be radiation, we shall construct scenarios of transition from tachyon driven inflation to 
radiation domination for the following three possible types of the decay rate V [cf. Eq. (8) 
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below] of the inflaton into radiation: (i) a constant, (ii) dependent on the tachyon, and 
(iii) explicitly dependent on time [11]. 

2.1. Background equations in the presence of interacting components 

Consider a (3 + 1) -dimensional, spatially flat, smooth, expanding, Friedmann universe 
described by the line element 

ds 2 = dt 2 -a 2 (t)dx 2 , (1) 

where t denotes the cosmic time and a(t) is the scale factor. Let p and p denote the total 
energy density and the total pressure of a system consisting of multiple components 
of fields and fluids that are driving the expansion. Then, the Einstein's equations 
corresponding to the above line-element lead to the following Friedmann equations for 
the scale factor a(t): 

H 2 = (^-) p and H = - (AttG) (p + p), (2) 



3 

where H — (a /a) is the Hubble parameter. Also, the conservation of the total energy 
of the system leads to the continuity equation 

p + 3H (p + p) = 0. (3) 

The total energy density and the total pressure of the system can be expressed as 
the sum of the energy density p a and the pressure p a of the individual components as 
follows: 



Pa 



and p = J2p a - (4) 



If the different components of fields and fluids do not interact, then, in addition to the 
total energy density, the energy density of the individual components will be conserved 
as well. Hence, in such a situation, the energy density p a of each component will 
individually satisfy the continuity equation 

p a +3H (p a +pj = 0. (5) 

On the other hand, when the different components interact, the continuity equation for 
the individual components can be expressed as (see, for instance, Refs. [11, 28, 29]) 

p a +3H ( Pa+Pa ) = Q a , (6) 

where Q a denotes the rate at which energy density is transferred to the component a 
from the other components. The conservation of the energy of the complete system then 
leads to the following constraint on the total rate of transfer of the energy densities: 

E Q« = °- ( 7 ) 
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2.2. The case of the tachyon and a perfect fluid 

Let us now consider the system of a tachyon (which we shall denote as T) that is 
interacting with a perfect fluid (referred to, hereafter, as F), so that a = (T,F). We 
shall assume that the rate at which energy density is transferred to the perfect fluid, 
i.e. Q F , is given by [24] 

Q F =Tf 2 p T , (8) 

where p T is the energy density of the tachyon. Such a transfer of energy is assumed to 
describe — albeit, in a course grained fashion — the perturbative decay of the tachyon into 
particles that constitute the perfect fluid. The quantity T represents the corresponding 
decay rate and, as we shall discuss below, it can either be a constant, depend on the 
tachyon, or depend explicitly on time [11]. The last possibility effectively describes a 
situation wherein the inflaton first decays into particles corresponding to an intermediate 
field before further decaying into particles constituting the perfect fluid. 

For the above Q F , the continuity equation (6) corresponding to the perfect fluid is 
given by 

p F +3H (l+w F ) p F =TT 2 p T , (9) 

where w F = (p F /p F ) is the equation of state parameter describing the perfect fluid, 
which we shall assume to be a constant. Also, it is evident from Eq. (7) that 

Q T = -Q F = -Tt 2 p T . (10) 

Therefore, the continuity equation for the tachyon energy density p T reduces to 

p T + 3H (p T +p T ) = -TT 2 p T , (11) 

where p T is the pressure associated with the tachyon. Given a potential V(T) describing 
the tachyon, the corresponding energy density p T and pressure p T are given by (see, for 
instance, Refs. [4, 21, 22]) 

p T = ( -S22= ) and p T = -V(T) Vl-fi. (12) 

On substituting these two expressions in the continuity equation (11), we arrive at the 
following equation of motion for the tachyon T [24]: 

f \ ■ (V \ 

+ 3HT + TT+ (-f) =0, (13) 



l-T 2 \V 

where V T = (dV/dT). 

Before we proceed, the following clarification on the choice of Q F in Eq. (8) is in 
order at this stage of the discussion. Recall that, in a situation wherein the canonical 
scalar field, say 0, drives inflation, it is common to introduce a (r<j>) term in the field 
equation to describe the perturbative decay of the inflaton [see Refs. [1, 2, 15, 16, 17], 
also see Eq. (A. 4)]. Motivated by the canonical case, the choice of Q F in Eq. (8) has 
been specifically made so as to lead to the (TT) term in the equation of motion (13) 
for the tachyon [24]. Needless to add, the choice (8) is but one of the many possibilities 
that can help in achieving reheating at the end of inflation. 
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2. 3. Tachyonic inflationary models, different possible V and reheating 

In this sub-section, assuming the perfect fluid to be radiation with w F = (1/3), we shall 
construct specific scenarios of transition from tachyon driven inflation to an epoch of 
radiation domination. 

We shall consider two different types of potentials in order to describe the tachyon. 

• The first potential we shall consider is given by 

a potential that is well motivated from the string theory perspective [19]. 

• Our second choice will be the following phenomenologically motivated power law 
potential that has been considered earlier in the literature (see, for example, 
Ref. [21]): 

UT) = (— -A-J) . (15) 



X 1 + (T/T )\ 

In order to achieve the necessary amount of inflation and the correct amplitude for the 
scalar perturbations, suitable values for the two parameters A and T that describe the 
above potentials can be arrived at as follows. Firstly, one finds that, in these potentials, 
inflation typically occurs around T ~ T corresponding to an energy scale of about A 1 / 4 . 
Secondly, it turns out that, the quantity (AT 2 /Mp ) has to be much larger than unity 
for the potential slow roll parameters to be small and thereby ensure that, at least, 
60 e-folds of inflation takes place. The COBE normalization condition for the scalar 
perturbations essentially constrains the value of A. We shall then choose a sufficiently 
large value of (AT^/MpJ, thereby determining the value of T . 

In the absence of the fluid, we find that the two potentials V\(T) and V^(T) above 
allow about 60 e-folds of slow roll inflation and lead to the correct COBE amplitude 
for the choice of the parameters and initial conditions listed in Table 1. Also, in such a 



Potential 


A 


T 


(Ti/T ) 


T 


Vi{T) 


10~ 6 


10 5 


5.81 


0.1 


V2(T) 


10 -5.15 


3.7 x 10 4 


4.56 


0.1 



Table 1. Achieving inflation: The values for the parameters of the tachyon potential 
and the initial conditions for the inflaton (the initial value Tj of the tachyon and 
its initial velocity Tj) that we work with. These choices for the parameters and 
initial conditions lead to the required 60 e-folds of slow roll inflation and the observed 
amplitude for the scalar perturbations. 



situation, it is known that, at the end of inflation, the tachyon leads to an epoch of dust 
like behavior [20, 21], thereby making a radiation dominated epoch difficult to achieve. 
However, as we shall illustrate below, when the tachyon is interacting with the fluid, 
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we can ensure that a transition from inflation to the radiation dominated regime occurs 
with a suitable choice of the amplitude for the decay rate T. 

As we mentioned before, we shall consider the following three possible choices for 
the decay rate T: (i) I\ = A = constant, (ii) T 2 = T(T), and (hi) T 3 = T(t). We shall 
work with the following specific forms of T(T) and T(t): 

T 2 = T(T) = A (l - B tanh [(T - C) /£>]) , (16) 

T 3 = r(f) = A(l-B tanh [(* - C) /£>]) , (17) 

where A, B, C and T> are constants that we shall choose suitably to achieve the desired 
evolution. (It turns out that the above forms of T(T) and T(t) prove to be convenient 
in the numerical calculations, helping us illustrate the required behavior. We shall 
comment further on these choices in the concluding section.) We find that a transition 
from inflation to radiation domination can be achieved provided the amplitude of T 
(viz. A) is chosen to be less than the Hubble parameter H (which is almost a constant) 
during the epoch of slow roll inflation. In Table 2, we have listed the values of the initial 
energy density of radiation (viz. p l ) and the various parameters describing the decay 
rates (for the tachyon, we have used the same values listed in Table 1) that lead to 
about 60 e-folds of slow roll inflation, the correct scalar amplitude, while also reheating 
the universe within 2-3 e-folds after the end of inflation. In the left column of Figure 1, 



Potential 


r 


p\ 


A 


B 


C 


V 


Vi(T) 




1Q -40 


5A 








r 2 


1Q -40 


5A 


io- 2 


2T 


T 


r 3 


10 -40 


50 A 


-0.98 


3.09 x 10 6 


10 5 






1Q -40 


1.43 x 10" 2 A 








r 2 


1Q -40 


10" 2 A 


io- 2 


2T 


T 


r 3 


1Q -40 


lO^A 


-0.98 


3.29 x 10 6 


10 5 



Table 2. Achieving inflation and reheating: The values of the initial energy density 
of radiation p* and the parameters describing the decay rates that we work with when 
the interaction of the tachyon with radiation is taken into account. For the tachyon, 
we use the same values listed in the previous table. These parameters and initial 
conditions lead to about 60 e-folds of slow roll inflation, the required amplitude for the 
perturbations, and also reheat the universe within a couple of e-folds after inflation. 

we have plotted the fractional contributions (with respect to the total energy density) 
of the energy densities of the tachyon and radiation, i.e. Q T = (p T /p) and f2 7 = (p 7 /p), 
as a function of the number of e-folds N for the potential V\(T). It is evident from the 
figure that the various choices for the decay rate V ensure the completion of reheating 
(i.e. n ~ 1) within a couple of e-folds after the end of inflation. In the right column 
of the figure, we have plotted the corresponding effective equation of state parameter, 
viz. w = (p/p), and also the first Hubble slow roll parameter, viz. e = —(H/H 2 ), of 
the entire system, for the same tachyon potential. These plots, while confirming that 
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N N 

Figure 1. In the left column, the evolution of the quantities Q T (in blue) and f2 T 
(in red) has been plotted as a function of the number of e-folds N for the the tachyon 
potential V\(T) and the three cases of V, viz. Ti = constant (on top), T2 = T(T) 
(in the middle) and T3 — T(t) (at the bottom). Similarly, in the right column, the 
evolution of the first slow roll parameter e (in blue) and the equation of state of the 
entire system w (in red) has been plotted as a function of the number of e-folds for 
the same potential and the different cases of T, as indicated above. All these plots 
correspond to the values for the various parameters and the initial conditions that we 
have listed in Tables 1 and 2. The figures clearly illustrate the transfer of energy from 
the inflaton to radiation. (Note that e = 2 during the radiation dominated epoch.) We 
have chosen the various parameters in such a fashion that there is a sufficiently rapid 
transition from inflation to radiation domination without any intermediate regime, as 
suggested by the evolution of e and w. We find that a very similar behavior occurs for 
the potential V 2 (T). 
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inflation has indeed ended and reheating has been realized, also indicate the nature of 
the composite matter during the transition. We should add that a very similar behavior 
occurs for the potential V^T). 



3. Equations of motion governing the scalar perturbations for the system 

In this section, using the first order Einstein's equations and the equations describing 
the conservation of the perturbed energy density of the tachyon and the perfect fluid, 
we shall arrive at the coupled, second order (in time) differential equations governing 
the scalar perturbations for the system consisting of the tachyon and the fluid. For 
convenience, we shall work in the uniform curvature gauge (UCG), a gauge that is also 
referred to, often, as the spatially flat gauge. 



3.1. First order Einstein's equations 

If we take into account the scalar perturbations to the spatially flat background 
metric (1), then, in the UCG, the Friedmann line-element is given by [1, 2] 

ds 2 = (1 + 2 A) dt 2 - 2 a (<%£) dxi dt - a 2 (t) dx 2 , (18) 

where A and B are functions that describe the two degrees of freedom associated with 
the perturbations. Being a scalar field, the tachyon, evidently, does not possess any 
anisotropic stress at the linear order in the perturbations. We shall assume that the 
perfect fluid does not contain any anisotropic stress either. In such a case, at the first 
order in the perturbations, the Einstein equations in the UCG can be written as [1, 2] 

-(3H 2 )A- (tpj V 2 B= (4nG)5p, (19) 

H(d t A) =(4ttG) (diSq), (20) 

HA+ (2H + 3# 2 ) A =(4irG)5p, (21) 

where the quantities Sp, (diSq), and 5p denote the perturbations in the total energy 
density, the total momentum flux, and the total pressure of the complete system, 
respectively. The absence of the anisotropic stress leads to the following additional 
relation between the two functions A and B describing the scalar perturbations: 

A + a (b + 2HB > ) =0. (22) 

For the system of our interest, viz. that of the tachyon and a perfect fluid described by 
a constant equation of state parameter w F , the quantities 5p, 5q and 5p are given by 

5 P = (5p T + 5p F ) = \ '^——^ ) + ; yT m (ST -At)+ Sp F , (23) 
Sq = (5q T + 6q F ) = I — I - ^, (24) 
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5p = (Sp T + 5p F ) = - (y T ST Vl - T 2 ) + ^-J=J-^ (ST -AT^+w F Sp F , (25) 

where ST denotes the perturbation in the tachyon and ip F is proportional to the potential 
that determines the three (i.e. the spatial) velocity of the perfect fluid:):. 



3.2. Equations of motion describing the perturbed matter components 

In the UCG, at the first order in the perturbations, the equation describing the 
conservation of the energy density of the component a is given by [2, 27, 28, 29] 

Sp a + 3 H (S Pa + S P J- (J^j V 2 [Sq a + ( Pa+ p a )aB]-Q a A- SQ a = 0, (26) 

where the quantities Sp a , Sq a and Sp a denote the perturbations at the linear order in the 
energy density, the momentum flux, and the pressure of the particular component a, 
while Q a and SQ a indicate the rate at which the energy is transferred to the component a 
and the perturbation in the rate, respectively. The equation of motion that governs the 
perturbation ST in the tachyon can now be obtained upon using the expressions (23)- 
(25) for Sp T , Sq T and Sp T in above equation for the conservation of the perturbed energy 
density. We find that the equation of motion describing the perturbation ST is given by 




-IT 



ST 



V 



ST 



1 — T 2 



6HT +2 -f +TT 1 +T 



1 — T 2 



- - (V 2 B) +TST = 0, 



(27) 



where ST denotes the first order perturbation in the decay rate. Upon using Eq. (26), 
we find that the equation describing the conservation of the perturbed energy density 
of the fluid Sp F can be written as 



Sp F +3H (1 



s Pf + 3 



^ YVT(2-f 2 ^ 



(l - T 2 ) 

(1 +w F )p 1 



3/2 



5T- 



I 



TV T T 2 
\/l-T 2 



ST + 



TVT 2 



(l - T 2 ) 



3/2 



A 



V 2 B- 



VT 2 



V 7 ! - t 2 



ST = 0. 



(28) 



t Usually, for a perfect fluid, the quantity Sq F is written as [{p F +p F )x F }> where x F is the potential 
whose spatial gradient describes the three velocity of the fluid (see, for instance, Rcfs. [1, 2]). For 
convenience, we have instead defined the entire quantity [(p F +p F ) x F \ as ~' i Pf [16]- As we shall point 
out later, we have ensured that suitable sub- Hubble initial conditions are imposed on ip F . 
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Also, on using the Einstein's equations (20) and (21), we obtain the following first order 
(in time) differential equation for the quantity ip F that describes the spatial velocity of 
the fluid: 

if> F + 3Hi/> F +w F 6p F + (l+w F )p F A+ I P=L ) 5T = 0. (29) 



1 — T 2 / 

We should emphasize here that the above equations take into account all the three 
possibilities of T (i.e. it can either be a constant, depend on the tachyon, or on time) 
that we had discussed earlier. 



3.3. The coupled perturbation variables and the initial conditions 



The variables describing the perturbations in the tachyon and the fluid that we shall 
eventually evolve numerically are [2, 16] 



Q T = 5T and Q f 



1>* 



(30) 



We have already obtained a second order differential equation (in time) describing 
the evolution of Q T [viz. Eq. (27)]. A similar equation for Q F can be obtained by 
differentiating Eq. (29) with respect to time, and upon using Eq. (28) that describes 
the conservation of the perturbed energy density of the fluid. Then, on utilizing 
the first order Einstein's equations (19), (20) and (22) to eliminate the scalar metric 
perturbations A and B, we can arrive at a set of coupled, second order (in time, again!) 
differential equations for the quantities Q T and Q F . 

For simplicity, let us first discuss the case wherein the decay rate T is a constant. 
We shall later indicate as to how certain coefficients in the differential equations change 
when T is assumed to be either dependent on the tachyon or on time. Upon suitably 
using the background equations (2), (9) and (13), we find that the differential equations 
satisfied by the perturbation variables Q T and Q F can be written as 



(31) 



where, as before, (a,/3) = (T,F). The coefficients J r a/3 and Q af} appearing in the above 
equations are given by 



F TT = (3H + T) (l-3T 2 )-2T (^j , 



*/ 1-7-1 



** IT 1 
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(39) 



3HT 
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TVT 3 



As we had pointed out, the coefficients J r afj and Q afj above have been derived 
assuming that the decay rate T is a constant. For the case wherein the decay rate 
is a function of the tachyon, the coefficients Q TT and Q FT are modified to 

Q TT — > G TT + (l — T 2 j f T 2 
VT 2 



Oft + 



(40) 
(41) 



\/l —f 2 ) \y/(l+W F ) p 

(dT/dT). Similarly, in the case wherein the decay rate is an explicit 



where T T 

function of time, we find that these two coefficients are modified as follows: 




(42) 



(43) 



The other coefficients remain unaffected in both these cases. 

We should mention here that we have checked that we indeed recover the standard 
equations in the different limiting cases from the coupled equations (31) for Q T and Q F . 
For instance, when the fluid is ignored and the decay rate is assumed to vanish, the 
equation for the variable Q T reduces to the equation that has been derived earlier 
in the literature for describing the perturbation in the tachyon (see, for instance, 
Refs. [4, 21, 25]). Moreover, when the tachyon is assumed to vanish and its coupling to 
the fluid is also ignored, we find that we arrive at the well known equation governing 
the fluid perturbation, as required [1, 2]. 



4. Effects of reheating on the scalar power spectrum 



In this section, we shall discuss the effects of reheating on the large scale curvature 
perturbations and the scalar power spectrum. We shall also investigate the behavior 



Reheating in tachyonic inflation and the effects on the large scale perturbations 



15 



of the non-adiabatic pressure perturbations at large scales, when reheating is achieved 
through the three different choices for the decay rate. 



4-1. Evolution of the curvature perturbations and the scalar power spectrum 

We shall solve the coupled equations (31) for the two variables Q T and Q F numerically. 
Imposing the standard Bunch-Davies initial conditions, we shall evolve the perturbations 
from the sub-Hubble to the super- Hubble scales, and analyse the effects of the transition 
from inflation to the radiation dominated epoch on the evolution of the large scale 
perturbations. We shall compute the scalar power spectrum at the end of inflation and 
soon after reheating is complete. Let us begin by quickly summarizing the essential 
quantities of interest and the initial conditions that we shall be imposing. 

Let ft T and 1Z F denote the curvature perturbations associated with the tachyon 
and the fluid, respectively. We find that these curvature perturbations are related to 
the variables Q T and Q F as follows [2, 4, 21, 25]: 



ft. 



ft. 




H 



Q, 



(44) 



(45) 



.Pf+PfJ " \y/(l+W F ) p F 

The total curvature perturbation of the system, say, ft, can then be expressed as a 
weighted sum of the individual curvature perturbations in the following fashion (see, for 
instance, Refs. [28, 29]): 

'Pa +Po 



ft 



E 



P + P 



ft. 



(46) 



Therefore, for our system of interest, the total curvature perturbation is given by 

-1 



ft = H 



VT 2 



Vi-f 2 
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VT 



Q T - y/[{l+W F ) p F ] Q f 



(47) 



1 -T\ 

The scalar power spectrum is then defined in terms of the Fourier mode TZk of the total 
curvature perturbation as 

k 3 



2tt 5 



Iftk 



(48) 



The Mukhanov- Sasaki variables, say, v T and v F , associated with the tachyon and 
the fluid are related to the corresponding curvature perturbations as follows: 

v T = (ft T z T ) and v F = (ft F z F ) , (49) 

where the quantities z T and z F are given by [1, 2, 4, 21] 



8ttG 



1/2 



aT 



and 



nl/2 



(50) 
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As is usually done, we shall impose the following initial conditions on the Fourier modes 
of the perturbation variables v T and v F : 

V2 /„-\ .,,.,1/2 



2uj k I \a J \ 2 / 



and * = - i (51) 



where 



((A;c s ) 2 -a 2 [(^) + (#*/*)]), (52) 

with Cg = (1 — T 2 ) in the case of the tachyon and c 2 s = w F in the case of the fluid. 
Also, we shall impose these conditions when the modes are well inside the Hubble 
radius during the inflationary epoch. We shall transform these initial conditions to the 
corresponding conditions on the variables Q T and Q F when numerically integrating the 
coupled equations (31). 

We shall focus on modes that correspond to cosmological scales today. We shall 
choose our parameters in such a fashion that modes spread over four orders of magnitude 
(say, 10 -4 < k < 1) leave the Hubble radius during the early stages of inflation (they 
actually leave during 5 < N < 15). As we have discussed earlier, for our choice of 
the parameters and initial conditions on the background variables, reheating is achieved 
after about 65 odd e-folds (also see Figure 1). Moreover, the transition from inflation to 
radiation domination occurs within about 2-3 e-folds after inflation has terminated. We 
evaluate the scalar power spectrum towards the end of inflation as well as soon after the 
transition to the radiation dominated epoch is complete. In the left column of Figure 2, 
we have plotted the evolution of the amplitudes of the individual (viz. \7Z T \ and \7l \ 
corresponding to the tachyon and radiation, respectively) as well as the total curvature 
perturbation (i.e. \1Z\) for a typical cosmological scale (we have chosen k — 0.1) across 
the transition from inflation to radiation domination. And, in the right column of the 
figure, we have plotted the scalar power spectrum (constructed out of the total curvature 
perturbation) evaluated before and after the transition. It is evident from the plots that 
the amplitude of the total curvature perturbation remains unaffected during reheating 
when the decay rate T is either a constant or is a function of the tachyon. However, 
when the decay rate T is an explicit function of time, for our choice of T(t), we find that 
the amplitude of the total curvature perturbation is enhanced during reheating. These 
results confirm similar conclusions that have been arrived at earlier based on the first 
order (in time) super-Hubble equations [11]. As we shall illustrate in the next section, 
this enhancement can be attributed to the growth of the tachyon-fluid non-adiabatic 
pressure perturbations during the transition. 



4-2. Entropy perturbations as the source of non-trivial evolution 

In this section, we shall discuss the evolution of the entropic (i.e. the non-adiabatic 
pressure) perturbations across the transition from inflation to the radiation dominated 
epoch. We shall show that it is the growth of the isocurvature perturbations (i.e. 
relative non-adiabatic pressure perturbation between the tachyon and the fluid) that is 
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Figure 2. In the left column, the evolution of the amplitudes of the curvature 
perturbations associated with the tachyon (\1Z T \, in blue) and radiation (|72. |, in red) 
as well as the weighted sum of the two (\1Z\, in green) has been plotted as a function 
of the number of e-folds N for the three possible cases of T, viz. Ti = constant (on 
top), T 2 = T(T) (in the middle) and T 3 = T(t) (at the bottom). In these plots, we 
have chosen a typical cosmological scale with the wavenumber k = 0.1. Clearly, the 
amplitude of the total curvature perturbation remains constant across the transition in 
the first and the second cases, while in the third case, it is enhanced during reheating. 
In the right column, the scalar power spectrum has been plotted towards the end of 
inflation (in blue) and soon after complete reheating has been achieved (in red) for 
all the three cases of T, as mentioned above. These plots unambiguously illustrate 
that the amplitude of the nearly scale invariant power spectrum is modified only in 
the case wherein the decay rate is explicitly time dependent. Note that all these plots 
correspond to the potential V\ (T) with values for the various parameters and the initial 
conditions that are listed in Tables 1 and 2. We find that a very similar behavior occurs 
for the potential V 2 (T). 
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responsible for the change in the amplitude of the large scale curvature perturbations 
during reheating. 

Let us recall the key quantities and equations. In the UCG, the gauge invariant 
Bardeen potential <E> is given by 

$ = A + a (b + HB?) . (53) 

Upon using the first order Einstein equations (19), (21) and (22), it can be shown that 
the Bardeen potential satisfies the differential equation [1, 2] 



$+ (4 + 3c 2 ) H&+ 2H + 3H 2 (l + c 



-A V 2 $ = (AttG) 5p NA , (54) 



where c 2 A = (p/p) denotes the adiabatic speed of sound in the composite system, and 
the quantity 5p NA is the non-adiabatic pressure perturbation defined as 

Using the Einstein equations (20) and (22), the total curvature perturbation 1Z of the 
system can be written in terms of the Bardeen potential $ as follows [1, 2]: 



K = *+ ^ • (56) 




On substituting this expression for the curvature perturbation in the equation (54) that 
governs the evolution of the Bardeen potential <E> and making use of the background 
equations, one can arrive at the well known equation 

h \ r / c 2 \ 

' ' (47rG) 5p NA + ( ) V 2 $ 



n = - 



H 



0/ 



(57) 



This equation clearly suggests that, at super-Hubble scales, any change in the amplitude 
of the curvature perturbations has to be due to the non-trivial evolution of the non- 
adiabatic pressure perturbation 5p KA . 

The perfect fluid, by definition, does not possess any intrinsic entropy perturbation. 
Therefore, the non-adiabatic pressure perturbation 5p NA of our composite system 
depends on the intrinsic non-adiabatic pressure perturbation of the tachyon, say, Sp^ A , 
and the relative non-adiabatic pressure perturbation between the tachyon and the 
perfect fluid (i.e. the isocurvature perturbation), which we shall denote as 5p^ A . Hence, 
we can write 

5p NA = (5p N T A + 6p" T A ) , (58) 
where the quantity Sp^, A is defined as (see, for instance, Ref. [4]) 

5p» A = {5p T - c 2 5p T ) (59) 
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with c 2 = (p T /p T ). The relative non-adiabatic pressure perturbation Sp^ p is then given 
by (see, for instance, Refs. [11, 30, 31]) 

' Pt +Pt 



5p^ = (3H + T) 



p + p 



[(p F+ p F )-(T/3H) (p T +p T )} 



5 Pl 



5p i 



(60) 



x (w F -c 2 ) 

V\p T ) \Pf 

Upon using the expressions (23) and (25), and the Einstein's equations (19), (21) 
and (22), we find that we can write the quantities 5p T , 5p T and 5p F in terms of the 
variables Q T and Q F as follows: 
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(63) 



These three relations allow us to express the non-adiabatic pressure perturbations 5p^ A 
and 5p£A in terms of Q T , Q F , and their time derivatives. 

In Figure 3, we have plotted the amplitudes of the Bardeen potential (i.e. |$|) and 
the non-adiabatic pressure perturbations (viz. |5p^ A | and |5p£ A |) as a function of the 
number of e-folds. For convenience in numerical computation, we have chosen a very 
small scale mode (corresponding to k — 10 16 ) that leaves the Hubble radius around 
58 e-folds or so during inflation§. The following points are evident from the three 

§ The amplitude of the non-adiabatic pressure perturbations corresponding to the cosmological scales 
(i.e. 10~ 4 < k < 1) prove to be very small when they are well outside the Hubble radius. As a result, 
evolving these reliably until the transition to the radiation dominated epoch requires considerable 
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Figure 3. The evolution of the amplitudes of the Bardeen potential |$| (in green), 
the intrinsic entropy perturbation |<5p~ A | (in blue) and the isocurvature perturbation 
\Sp^\ (in red) has been plotted as a function of the number of e- folds N for the three 
different cases of T, as in the previous two figures. For computational reasons, we 
have chosen a very small scale mode with wavenumber k — 10 16 that exits the Hubble 
radius during inflation around 58 e-folds or so. Also, we have normalized the values 
of all the three quantities to be unity at Hubble exit. Clearly, while the Bardeen 
potential remains nearly a constant across the transition, both the intrinsic entropy 
and the isocurvature perturbations decay in almost the same fashion before as well 
as after the transition. The dashed lines in black and cyan indicate the e _ ' 2JV ' and 
e -(80JV) behavior, respectively. Interestingly, the non-adiabatic pressure perturbations 
die down extraordinarily rapidly after reheating is complete. The third plot illustrates 
the growth in the isocurvature perturbation for a short duration before the transition. 
As we have discussed in the text, it is this growth that is evidently responsible for the 
amplification of the curvature perturbations when the decay rate is explicitly dependent 
on time. Note that all these plots correspond to the potential V\(T) with values for 
the various parameters and the initial conditions that we have used in the last two 
figures. As in the earlier figures, we find that a very similar behavior occurs for the 
potential V 2 (T). 
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plots. To begin with, barring a very small change in its amplitude — which is expected 
during the transition from inflation to the radiation dominated epoch — the Bardeen 
potential <£> remains a constant at super-Hubble scales. Secondly, it is clear from the 
plots that the intrinsic entropy perturbation associated with the tachyon decays as e~( 2N ^ 
at super-Hubble scales during inflation, a result that is well known in the literature 
(see, for instance, Refs. [3, 4]). Thirdly, and interestingly, we find that the isocurvature 
perturbation behaves in almost the same fashion as the intrinsic entropy perturbation 
both before and after the transition, a feature that has been noticed earlier [31]. While 
they both behave as e~^ 2N ^ at super-Hubble scales before the transition, they die down 
extremely rapidly (roughly as e~^ 80Ar ^!) during the radiation dominated epoch. Lastly, it 
is also evident from the third plot that, when the decay rate is a function of time, while 
the intrinsic entropy perturbation continues to decay, the isocurvature perturbation 
exhibits a brief period of growth close to the transition before decaying. Obviously, it 
is this growth of the isocurvature perturbations that is responsible for the amplification 
of the curvature perturbations during the transition. 

5. Summary and discussion 

In this work, we have studied the evolution of perturbations in an interacting system 
consisting of a tachyon and radiation. Treating the tachyon as an inflaton, we have 
investigated the effects of reheating — i.e. the perturbative transfer of energy from the 
tachyon to radiation — on the large scale curvature perturbations. We have shown that 
the transition does not alter the amplitude of the total curvature perturbation of the 
system when the rate describing the decay of the inflaton into radiation is either a 
constant or a function of the tachyon. However, we find that the total curvature 
perturbation is modified when the decay rate associated with the tachyon turns out 
to be explicitly dependent on time. We find that, for the specific functional form of the 
decay rate that we work with, the overall amplitude of the power spectrum is enhanced 
(by a factor of about 25) during the transition. We have also been able to explicitly 
illustrate that it is the growth of the isocurvature perturbations close to the transition 
is responsible for the amplification of the curvature perturbations. 

It may be considered an overkill to integrate the second order differential equations 
rather than work with the much simpler first order equations [11, 27, 28]. Our 
motivations were threefold. Importantly, our effort allowed us to cross-check certain 
previous results [11]. Also, our aim is to later study the effects of deviations from slow 
roll inflation as well as the effects of the modified dynamics close to the transition to the 
radiation dominated epoch on the small scale perturbations. As we had mentioned in the 

numerical accuracy, which is difficult to achieve. For this reason, to illustrate the evolution of the 
non-adiabatic pressure perturbations, we choose to work with a very small scale mode that leaves the 
Hubble radius close to the end of inflation. Since the inflationary epoch is of the slow roll type for 
all the e-folds, the non-adiabatic pressure perturbations associated with the cosmological scales can be 
expected to behave in a fashion similar to that associated with such a small scale mode. 
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introduction, these effects can have interesting implications on the number of primordial 
black holes that are formed towards the end of inflation [7]. Moreover, in possibilities 
such as the warm inflationary scenario, the decay rate, in addition to depending on the 
inflaton, can also depend on the temperature (see, for example, Refs. [25, 26, 32]). In 
such a case, the transfer of the energy from the inflaton to radiation may affect the 
amplitude and possibly the spectral index as well [17]. 

We would like to conclude by commenting on our choices for T(T) and V(t). As 
we had pointed our earlier, our choices were motivated by convenience in numerical 
evolution rather than physics. Obviously, better motivated functional forms of these 
quantities need to be investigated. In particular, it will be interesting to understand 
whether there exist T(t) that suppress the amplitude of the curvature perturbations 
during the transition rather than enhance them. Such a possibility is worth exploring 
as it can lead to an increase in the tensor-to-scalar ratio post inflation (see Ref. [33]; 
however, see also Ref. [34]). We are currently investigating these issues. 
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Appendix A. The case of the canonical scalar field and a perfect fluid 

In this appendix, we shall illustrate the corresponding effects for the case of the 
canonical scalar field model that was considered recently in the literature [16]. We 
shall rapidly summarize the essential equations describing the background evolution 
and the perturbations, and present the results. 

Appendix A. 1. Background equations 

In the case of a canonical scalar field, say, 0, that is interacting with a perfect fluid, if 
one assumes that Q F = (r<^ 2 ), then the equation (6) that describes the conservation of 
the energy density of the fluid is given by 



p F + 3H (l + w F ) p F = r> 2 . 
Then, Eq. (7) implies that = -Q F = -(T 2 ) and, 




governing the energy density p of the scalar field reduces to 



(A.2) 
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where p^ denotes the pressure of the field. The energy density and pressure associated 
with the canonical scalar field are given by 



Upon using these expressions in the continuity equation (A. 2), we arrive at the following 
equation of motion governing the scalar field [15, 16, 17]: 

+ 3#0 + T + 1^ = 0, (A.4) 

where = (dV/d(f>). 

Appendix A. 2. Equations governing the scalar perturbations 

Upon using the expressions for the perturbed energy density, momentum flux and 
pressure of the scalar field in Eq. (26), we can arrive at the equation of motion for 
the perturbation in the scalar field, say, 5(f), exactly as we did for the perturbation in 
the tachyon. In the UCG, we find that 5(f> satisfies the differential equation 

6<j> + (3H + r)6<j>+ v; 

-<j>A + (2V; + I>J A- I (V 2 B)+(f)5T = 0, (A.5) 

where = (d 2 V/d(f> 2 ) and, as before, 5T denotes the first order perturbation in 
the decay rate. Again, on using Eq. (26), we find that the equation describing the 
conservation of the perturbed energy density of the fluid 5p F can be written as 




5p F + 3H{l + w F )5p F +(\)V 2 ^ f 



-2F(f)5(f) + r ( f) z A- 



[l + w F )p f 



(V 2 B) - 2 5T = 0. (A.6) 



Also, from the first order Einstein's equations (20) and (21), we can arrive at the 
following differential equation for the quantity if) F that describes the spatial velocity 
of the fluid: 

i) F +3Hip F +w F 5p F + (l+w F ) p F A + T(p5(i) = 0. (A.7) 

The perturbation variables that we shall evolve numerically are = 5(f) and the 
quantity Q F that we had introduced earlier [cf. Eq. (30)]. As in the case of the tachyon 
and the perfect fluid, using the background equations, the first order Einstein's equations 
for the system and the energy conservation equations for the scalar field and the fluid, 
we can arrive at the coupled second order differential equations (31) for <2 and Q F , 
with (a, (3) = ((f>, F) [16]. For the case wherein the decay rate T is a constant, the 
coefficients J- ' ajj and Q afj are now given by 

F„ = (3H + r), (A.8) 
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For the case wherein the decay rate is a function of the scalar field, the coefficients 



and Q F<t> are modified to 



V( 1+W f) Pf 



(A.16) 
(A.17) 
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where = (dT/d0). In a similar fashion, when the decay rate is an explicit function 
of time, we find that the above two coefficients are modified as follows: 



The other coefficients remain unchanged. 
Appendix A. 3. Results 

In Figure Al, we have illustrated the transition from inflation to the radiation dominated 
epoch by plotting the dimensionless energy density parameters and f2 7 for the three 
possible types of decay We have also plotted the evolution of the first slow roll parameter 
e and the equation of state parameter w for the entire system. Note that we have 
considered the quadratic potential V(4>) = (m 2 4> 2 /2), and have chosen suitable values 
for the parameters and initial conditions to achieve the desired behavior. And, in 
Figure A2, as in the case of the tachyon discussed in the text, we have plotted evolution 
of the amplitude of the individual and the total curvature perturbations, as well as the 
scalar power spectrum before and after the transition. It is evident from these two 
figures that the conclusions we had arrived at for the tachyon apply equally well to the 
case involving the canonical scalar field as well. 
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Figure Al. In the left column, the evolution of the quantities SI, (in blue) and SI (in 
red) has been plotted as a function of the number of e-folds N for the three types of T, 
as discussed in the earlier figures. In the right column, we have plotted the evolution 
of the first slow roll parameter e (in blue) and the equation of state parameter w of 
the entire system. We have considered the popular quadratic potential to describe 
the scalar field, and we have worked with suitable parameters and initial condition to 
arrive at the required behavior. These plots clearly illustrate the transfer of energy 
from the inflaton to radiation. We should point out that the oscillations in the plots 
arise due the oscillations by the field at the bottom of the potential, a feature that 
does not arise in the case of the tachyon. 
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Figure A2. In the left column, the evolution of the amplitudes of the curvature 
perturbations associated with the scalar field (|7£J, in blue), radiation (\TZ |, in red) 
and the total curvature perturbation of the entire system (\TZ\, in green) has been 
plotted as a function of the number of e-folds N for the different types of the decay rate. 
For illustration, we have again chosen a typical cosmological mode with wavenumber 
k = 0.01. In the right column, as before, we have plotted the scalar power spectrum 
before (in blue) and after (in red) the transition. It is evident that, as in the case of 
the tachyon, the amplitude of the nearly scale invariant power spectrum is modified 
only in the case wherein the decay rate is explicitly time dependent. 
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